


gl 











m 





I 




SAFETY AND STABILITY 
IN 

CONCRETE BARREL SHELL 
ROOF STRUCTURES 

David Frederick Kelley 

«* *4 



A Thesis 

Presented to the Faculty of Princeton University 
in Candidacy for the Degree 
Master of Science in Engineering 



Recommended for Acceptance by the 
Department of Civil Engineering and Operations Research 



October 1991 



T253933 



I hereby declare that I am the sole author of this thesis. 

I authorize Princeton University to lend this thesis to other institutions or 
individuals for the purpose of scholarly research. 



I further authorize Princeton University to reproduce this thesis by photocopying 
or by other means, in total or in part, at the request of other institutions or 
individuals for the purpose of scholarly research. 













51 






Princeton University requires the signatures of all persons using or 
photocopying this thesis. Please sign below, and give address and date. 




I 



Acknowledgments 



I would like to express my heartfelt appreciation to my advisor, Professor David 
R Billington. In my initial conversation with him, we were discussing possible 
thesis topics, and he asked me why I had come to Princeton. I wasn’t quite sure 
how to answer the question at the time, but I am sure now. I came for the 
opportunity to work with quality faculty members like him. People who know 
that learning is more than what you get from a book, and teaching is more than 
standing in front of a blackboard. Thanks for the inspiration and patience. 
Professor! 

I am also indebted to my committee members. Professors Robert Mark and 
Jiann-Wen Ju, who provided equal inspiration through their courses, and who 
read this thesis and offered their suggestions. 

I would also like to thank the Civil Engineer Corps of the U.S. Navy, for allowing 
me this opportunity. It’s not just a job! 

A special thank you goes out to the many friends I have met during the year. I 
wish you continued success. Especially to John Matteo, who helped me in 
reviewing many projects by providing his keen insight and thoughtful words of 
wisdom. 

Additionally, I would like to thank Captain and Mrs. G, and Jane and Martin, who 
all introduced me to Princeton and gave me their never ending support. 

Finally, I would like to thank my wife, Teri, who stood by me all the way, and put 
up with the person who possessed my body for the last two weeks of writing this 
thesis. 



L 



Abstract 



The debate between Anton Tedesko and Charles S. Whitney which 
occurred from the 1930’s through the 1950’s typifies the confusion among 
designers in the United States regarding thin shell concrete roof design. Each 
man thought his method was correct and designed structures constructed in 
America during the first half of the twentieth century. By taking a closer look at 
their debate, we can gain some insight into their methods of design. To resolve 
the conflict, we then apply modern methods of analysis to analyze a hangar 
model Whitney had presented in his articles. A full span analysis is performed 
using the finite element computer program P-FRAME. In addition, we address 
concerns which were not incorporated into the original analysis. We employ the 
methods of Milo S. Ketchum and Robert S. Rowe to compute deflection 
moments for the structure. In addition, we use Ketchum and Rowe’s work as 
background for developing the Initial Deflection Method of computing buckling 
safety factors. To validate the procedure, we compute buckling safety factors 
for a variety of structures and compare them to classical formulations. 
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Chapter One 
Introduction 



Over the last century, designers such as Eugene Freyssinet, Robert 
Maillart, Pierre Luigi Nervi, Felix Candela and Heinz Isler have brought the art 
of reinforced concrete design to its mature state. 

Isler’s thin shell concrete roofs cover many European structures. From 
tennis courts to gas stations, his shells provide practical, yet interesting, 
solutions to everyday roofing problems. 

In the United States, however, reinforced concrete design has not 
advanced as it has abroad. A reflection of this lack of progress can be seen in 
the content of basic design texts. During the 1950’s, in the well known text book 
by George Winter and Arthur H. Nilson, Design of Concrete Structures, an 
entire chapter was dedicated to arch and shell design. In the 1991 edition, the 

words arch and shell do not even appear in the index.'' Why has this type of 
design disappeared from our basic text books? 

A possible explanation is that a general confusion exists in America 
regarding shell behavior and because of this, key safety questions still remain 
unanswered. 

This confusion is clearly demonstrated in the debate between Charles S. 
Whitney (1902-1961) and Anton Tedesko (b. 1902) which took place in the 
1940’s and 1950’s concerning concrete barrel shell roof design. The design is 
a thin concrete barrel shaped shell with arch stiffeners spaced along the length. 
The debate focused on how to position the shell in relation to the arches. 

Tedesko’s opinion was that the shell should be positioned at the rib 
extremity. Whitney, on the other hand, believed that the shell should be located 
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at the mid-height of the rib. They debated in many engineering publications, but 
without resolution. 

Tedesko was an Austrian born engineer who had studied civil 
engineering at the Technological Institute in Vienna in the early 1920’s, and had 
learned thin-shell concrete roof design while working at the firm of Dyckerhoff 
and Widmann in Weisbaden. In 1932, because of prior work experience in the 
U.S., he was sent to work in America when his firm decided to expand its 

operations.^ Once there, he gained an affiliation with the Roberts and Schaefer 
Company in Chicago, and stayed on with them to do extensive thin-shell design 
work during the 1930’s and 40’s. Because of his effort in this capacity, he 

introduced thin-shell concrete roof structures to the United States.^ Examples 
of his work are the Ice Hockey Arena in Hershey, Pennsylvania which opened in 
1936 and the U.S. Navy Hangars at North Island in San Diego, California 
designed in 1941 . 

Whitney was an American born designer who wrote two major articles on 
arch design in the 1920’s. In the 1940’s, he developed a new desing method 
for barrel shell roofs that he called "... a novel feature. . . . which has important 

advantages.” ^ In his method, he placed a great emphasis on volume change 

moments, which are a function of the cross sectional moment of inertia. The 
moments of inertia would be less if the shell was located in the middle of the rib, 
thus, this was the better design. 

Whitney also designed structures which were constructed in the United 
States. An example of his work is the Field House at Syracuse University built 
in the 1950’s. 

We can see, then, that two very different methods of design existed 
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simultaneously. The debate between Whitney and Tedesko which started in the 
literature over 40 years ago, has not been resolved. In this thesis, we will first 
examine the different design methods and then utilize modern engineering 
tools to clarify the debate. 

We look at Whitney’s method by examining calculations which were 
prepared for a hangar model and presented in articles published during the 
1940’s and 1950’s. Tedesko’s rebuttal to one of the articles is also scrutinized 
to examine his ideas on the subject. 

In our modern analysis, we use the finite element method to analyze 
Whitney’s hangar model. Additionally, we employ the methods of Robert S. 
Rowe and Milo S. Ketchum to calculate stress amplification due to deflections in 
arches. Using these as a starting point, we develop a method of predicting 
buckling loads for arched structures. 
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Chapter Two 



Whitney’s “New Idea” 



Whitney began writing about concrete arch design as early as 1925 with 
his article, “Design of Symmetrical Concrete Arches” published in the American 
Society of Civil Engineers (ASCE) Transactions. Between 1932 and 1940 he 
was the Chairman of the American Concrete Institute (ACI) Committee 312 
which was attempting to establish standards for reinforced concrete arch 
design. His ideas appear in this committee’s reports published in 1932 and 
1940. 

It wasn’t until after this that he began writing about his method of 
concrete barrel shell roof design. He presented his ideas in three articles 
published between 1943 and 1955. 

The structure Whitney analyzes in developing his claims is a 220 foot 
clear span aircraft hangar model. The hangar roof is a parabolic barrel shell 
comprised of a 4 inch thick reinforced concrete slab with ribs spaced 20 feet 
center to center. The rib cross section varies from 18x32 inches at the crown to 
18 X 40 inches at each springing. The height of the rib center line above the 
supports at mid-span is 27.5 feet, and the roof is supported at each springing by 
concrete A-frames. Whitney assumes that the A-frames act to fix the ends of the 
barrel shell by restricting translational and rotational motion. Figure 1 is a 
longitudinal and transverse section of Whitney’s model. This structure was first 
presented in his 1944 article “Aircraft Hangars of Reinforced Concrete”.^ 

The article that we will focus on was published in the ACI Journal in June 
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1950 under the title, “Cost of Long Span Concrete Roof Shells". In this article, 

Whitney explains the advantages of his design method: 

“An important feature of this type of construction is the placing of the shell 
near the neutral axis of the hbs so that the ribs project about half above 
and half below the shell. The principal effects of this arrangement 
structurally are the elimination of edge stresses in the shell due to rib 
flexure and the reduction of the stiffness of the combined rib and shell 
with a corresponding reduction in volume change moments.” ® 



He also develops a chart which shows how his shell positioning reduces 
required rib size and thus, construction cost. 
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Figure 1 . Longitudinal and transverse sections of the 220 ft span hangar model 



His chart, shown in Figure 2, presents three different crown cross 
sections, each based on a 20 foot spacing of the arch ribs. The first cross 
section has the shell located at the mid-height of the rib, the second has the 
shell positioned at the top of a similar rib, and the third has the shell at the top of 
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an enlarged rib. The larger rib in the third cross section is necessary, according 
to Whitney, to provide the equivalent strength of the first cross section. 

To understand Whitney’s method of design, we will examine the numbers 
in his chart. Since Whitney did not publish detailed calculations, we must use 
information from his other articles and reports to estimate his results. 

As his first table entry, Whitney gives the moment of inertia at the crown. 
For the first section, it is computed using the 20 foot width and including 
reinforcing steel at the top and bottom of the rib. In the next two sections, he 
only uses a 14 foot width in his calculations. He explains that with the slab at 
the top of the rib, only 70% of the shell is effective. Whitney does not provide 
any background for this assumption in any of his published material, but we can 

verify it with a formula published in a 1990 textbook on concrete shell design.^ 
The effective overhang of one side of the shell, b^, is: 

be = 0.76 (rh)i 
r = L2- 

where: 8d with: L= Span 

d = Rise 

h = Shell thickness 

Using Whitney’s data for the hangar model: 

be = 6.5 ft 

This gives a total overhang of 13 ft. When the rib width is considered, the total 
effective width becomes 14.5 ft, which compares very well with Whitney’s 14 ft 
assumption. 

Using Whitney’s data at the crown for the hangar model, our computed 
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moments of inertia for the three cross sections are: 



Section 1 : 58,000 in'* 

Section 2: 117,200 in'* 
Section 3: 154,800 in'* 



These compare very favorably to Whitney’s numbers. Detailed calculations 
appear in Appendix A. 
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Moment due to 
live load 


1,710,000 in.lb. 


1,710,000 in.lb. 


1,710,000 in.lb. 


Moment due to 
volume change 


1,026,000 in.lb. 


2,042,000 in.lb. 


2,670,000 in.lb. 


Total moment 


2,736,000 in.lb. 


3,752,000 in.lb. 


4,380,000 in.lb. 


Maximum horiz- 
ontal thrust 


447,1201b. 


450,170 lb. 


508,300 lb. 



Figure 2. Whitney s Data Table 



The second entry is live load bending moment. Whitney uses 30 psf as 
the live load for all three cross sections. For a 20 foot width, this results in a 
distributed load of; 

30 psf X 20 ft = 600 ^ 
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In his 1925 article, Whitney derives formulas to compute moments in concrete 
arches for different loads with varying cross sections. From Figure 50 in the 

1925 article, the maximum positive live load moment at the crown is:® 

Ml = Ki p L2 

where: p = Uniformly distributed live load 

L= Span length in feet 

Kj = Factor from Figure 50 using entering 

arguments N and m 

_ Yo _ quarterpoint rise - midspan rise 
with: ^ midspan rise 



m 



Ic 

Is cos 0s 



with: c = crown values 

s = springing values 



Using Whitney’s data for the hangar model: 

m=0.59 and N = 0.25 
With these as entering arguments for table 50: 

Ki = 0.0049 

The maximum positive live load moment at the crown is: 

Ml = 142,300 ft-lb = 1,710,000 in-lb 

Whitney presents this value in the table for all three cases. Even though the 
cross sectional length for the second and third cases is less, he is assuming the 
effective cross section carries the full live load. 

To determine load positioning, Whitney uses influence lines. Figures 35- 
39 in his 1925 article give influence lines for values of “m” ranging from 0.15 to 
0.40 and “N” values ranging from 0.15 to 0.25. From these figures, we can 
extrapolate the load positioning necessary to produce maximum positive 
moment at the crown for m = 0.59 and N = 0.25, which is shown in figure 3. 
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The third value Whitney lists is the volume change moment. The values 
are based on “rib shortening and a temperature drop of 40 degrees F including 

the effect of shrinkage.”® 





Rib shortening results from compressive axial forces in the structure, and 
the 40 °F temperature drop accounts for the worst case combination of 
temperature change and concrete shrinkage. Obviously, shrinkage produces 
an outward thrust with corresponding negative moments at the supports. With a 
temperature drop, these two effects would add together and make sense when 
considered with the 30 psf live load Whitney is using, which is probably a snow 
load. 

In the “Aircraft Hangars of Reinforced Concrete” article, Whitney provides 
insight to his choice of a 40 ®F temperature drop. First of all, he states that 
because plastic flow reduces temperature and shrinkage effects, 60% of the 
maximum temperature range for a geographical area should be used for 
concrete arch calculations. In the article, he presents a table of temperature 
ranges for various locations in the United States. We will choose a value from 
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the table for a location in the south, since this is what Whitney implies he used. 
The temperature range of 95 °F for New Orleans, Louisiana is chosen. The 
design range is therefore: 

95 °F X 60% = 57 °F 



Assuming that it drops from the mean, our change in temperature will be half of 
the design range, or 27 °F. Whitney suggests that the stress caused by 
shrinkage in concrete can also be represented by a drop in temperature. Using 

his shrinkage value of 15°F results in a total temperature drop 42 °F. 

The volume change moments are computed by first calculating the 
thrusts due to rib shortening and temperature changes. The thrusts are 
multiplied by a function of the rise to compute moments. Using Whitney’s 1925 
article, Micalos derives formulas for hingeless (fixed) arches. If we assume a 
secant variation in the cross section, the formula for computing thrust due to rib 

shortening is:'*'' 



Hrs 



= 45. iL i 
4 Am 



where: = Crown moment of inertia 

h = Mid span rise 
Am= Mean rib area 
H = Dead and live load thrust 



Using Whitney’s data, the thrust due to rib shortening is: 

Hrs = 4,180 lb 



The thrust due to temperature changes is computed from:'*^ 



Ht = ^ aTE^^ 
4 h2 
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where: 



a = Coefficient of thermal expansion 
T = Temperature change 



Using Whitney’s data and his recommended coefficient of thermal expansion of 
5.5 xIO’® in/in/°F, the thrust due to temperature change is:^^ 

Ht = 5270 lb 

The crown moment due to the total volume change thrust is: 

M= i h (Hrs + Hi ) 

With Whitney’s data, the crown moment due to volume changes is: 

M = 1,039,000 in-lb 

This value compares well with Whitney’s value of 1,026,000 in-lb. Whitney does 
not explain his assuming a secant variation in cross section for this calculation. 

The table entries for the other two sections vary directly with the moments 
of inertia. Thus, the volume change moment at the crown is nearly doubled 
when the shell is shifted from the mid-height to the extremity, and it is increased 
even more for the larger rib. 

The next table entry, the total moment, is simply the addition of the two 
previously calculated values. Since the live load moment is the same for all 
three cross sections, the total moment varies only with the change in volume 
change moments. 

Finally, Whitney computes the maximum horizontal thrust caused by the 
dead and live loads. Whitney’s formula for computing thrust due to dead load 

is:^4 
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^ Wc (g- 1) 
r 



where: g = The ratio of springing to crown weight 

Wc = Crown weight in Ib/ft 
L, = One-half the span length 

k = cosh"^ (g) 

Using 150 Ib/ft^ as the weight of concrete, the dead load thrust is: 

Hd = 337,600 lb 

The live load thrust for a load that produces the maximum positive moment at 
the crown is:"*^ 






where: Cj = A constant computed from Figure 

51 using m and N 
r = Midspan rise 
p = Distributed live load 



For Whitney’s data: 

Hi = 62,800 lb 

The dead and live load total thrust is: 

H = 400,400 lb 

This value is 46,720 lb lower than the table value of 447,120 lb. Whitney must be 
considering the full span live load in his calculations. To compute the full span 
live load thrust, we will use the previously defined formula with a Q value for 

both maximum positive and negative crown moment live loads:"' ® 
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Hi = 132,0001b 



The total thrust would then be; 

H = 469,600 lb 

This is 22,500 lb higher than the table value. Even if we consider the negative 
thrust from the volume change calculations, the computed total thrust would still 
be 13,000 lb higher than the table value. Whitney does not give any explanation 
for this difference. 

In summary, through his table, Whitney shows the importance of volume 
change moments in concrete barrel shell roof design. According to his 
conclusions, placing the shell at the mid-height of the rib cuts the volume 
change moment in half and is the most efficient design. He also claims that with 
the shell located at the extremity, the rib must be increased by 50% to give the 
equivalent strength of a cross section with the shell at mid-height. He does not, 
however, present any calculations to support this claim. 
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Chapter Three 
Tedesko’s Design Ideas 



Whitney’s ideas did not go unchallenged. The discussions of his papers 
raised serious questions as to the validity of his claims. In presenting the 
alternate viewpoint, we will focus on two papers in particular. The first is the 
discussion of Whitney’s 1950 paper “Cost of Long-Span Concrete Roof Shells’’, 
and the second is the discussion of the 1940 report of ACI Committee 312, 

“Plain and Reinforced Arches” . Both discussions were written by Structural 
Engineers from the Roberts and Schaefer Company of New York City, who were 

under the direction of Anton Tedesko at the time."' ^ Therefore, we will consider 
the alternate viewpoint as Tedesko’s. 

Tedesko does not agree with Whitney’s claim regarding shell position. 
Tedesko believes, instead, that the shell should be positioned at the rib 
extremity. He develops an alternate chart which shows that the rib width can be 
decreased if the shell is moved to the top or bottom. He supports his arguement 
by looking at the stress distribution and by computing buckling safety factors. 

We will examine Tedesko’s chart. Figure 4, to gain more insight into his 
argument. Tedesko uses four crown cross-sections which he calls cases one 
through four. The first two cross sections are the same as presented in 
Whitney’s table. The third has the shell at the bottom of the 1 8 inch wide rib, 
and the fourth has the reduced rib with the shell at the top. 

The first five items are the same ones listed by Whitney. Tedesko carries 
the calculations a bit farther, however, by computing stress distributions, 
tension force taken by the reinforcing steel and a buckling safety factor. 
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Case 1 


Case 2 


It* 

Case 3 


Moment of inertia 
of crown section, 
in. 


58,000 


116300 


116,300 


69,700 


Moment due to 
live load, in. -lb 


1,710,000 


1,710,000 


1,710,000 


1,710,000 


Moment due to 
volume change, 
in. -lb 


1,026,000 


2,042,000 


2,042,000 


1,230,000 


Total moment, 
in. -lb 


2,736,000 


3,752,000 


3,752,000 


2,940,000 


Max. horizontal 
thrust, lb 


447,120 


450,000 


450,000 


408,000 


Concrete stress 
at top of arch 
psi 


-1,096 


-628 


-1,091 


-696 


Concrete stress 
at bottom or arch, 
psi 


413 


404 


-59 


663 


Total tension 
force in concrete 
to be taken by 
reinf., lb 


32,500 
(min. reinf.- 
2.9 sq. in.) 


44,500 
(min. reinf.- 
2.9 sq. in.) 


0 

(min. reinf.- 
2.9 sq. in.) 


46,000 
(min. reinf.- 
2.9 sq. in.) 


Relative buckling 
safety of arch — 
proportional to I 
(Dischinger’s 
method) 


8.7 


17.4 


17.4 


10.5 



Figure 4. TedesKo's data chart 



The moments of inertia for the first two crown sections are the same as 
Whitney’s. Tedesko apparently agrees with Whitney’s use of the reduced 
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effective width when the shell is moved to the rib extremities. The moment of 
inertia for the third cross section is the same as the second since the two are 
mirror images. Their only differences will be in the section modulus for the top 
and bottom fibers, but this will only affect the stress distribution. For the fourth 
cross section, Tedesko computes the moment of inertia using the full 20 foot 
width. Since the shell is positioned at the top of the rib, however, the reduced 
effective width should be 14 feet. Using the reduced effective width, we 

calculate the moment of inertia to be 65,300 in'^; 4,400in'^ lower than Tedesko’s 
table value. 

The moment due to live load is the same as Whitney’s and the same for 
all four cases. Tedesko is not challenging Whitney’s use of the full width loading 
on the reduced effective width, the amount of load used, or his load position. 

The moments due to volume changes vary directly with the crown 
moment of inertia as in Whitney’s table, and the values for the first three cross 
sections are the same as before. Tedesko is not questioning the method of 
computation. The moment for the fourth case is about 6% too high due to 
Tedesko’s full width moment of inertia. Using Whitney’s volume change 
moment as a starting point, we calculate the moment for the fourth cross section 
to be 1,147,000 in-lb. This is 83,000 in-lb lower than Tedesko’s value. 

As in Whitney’s table, the total moments are computed directly from the 
live load and volume change moments. The total for case four using an 
effective width of 14 ft is 2,857,000 in-ib. 

Tedesko’s table values for horizontal thrust for the first three cases are 
the same as Whitney’s. He apparently concurs with the calculations from 
Whitney’s 1925 article. The thrust for the fourth cross section is computed using 
Whitney’s formulas as well. We concur with this value as it is not affected by the 
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change in moment of inertia. 

The next data entry in Tedesko’s chart is the stress distribution across the 
crown cross section. This is computed using standard formulas for axial and 
bending stress and the appropriate section modulus. Tedesko assumes an 
uncracked cross section, therefore, an elastic analysis is implied. Since the 
horizontal thrusts are all compressive, and the bending moments at the crown 
are all positive, general formulas for computing stresses at the extreme fibers 
are: 



top: ' A St bottom: A Sb 

where: P = Horizontal Thrust 

A = Cross sectional Area 
M= Total Bending Moment 
S = Section Modulus (I /y) 

Using Whitney’s data for the cross section in case one: 

ft = -1060 psi and fb= + 449psi 

Tedesko’s values are 36 psi lower for both the top and bottom fibers for case 
one. To arrive at the stresses listed in the table for the first cross section, 
Tedesko is either using a total thrust value of 499,800 lb instead of the listed 
value of 447,120 lbs or using a reduced cross sectional area. He does not 
explain the change. The stress calculations for the other three cases also differ 
from what would be expected from Whitney’s data, but they do show the trend 
Tedesko wants to demonstrate. In case four, the stresses are almost equally 
distributed and even with the reduced rib, the stresses are still well within the 
strength of concrete. 

Next, Tedesko lists the total tension force to be taken by the reinforcing 
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steel. This is calculated assuming a linear stress distribution and assuming the 
concrete does not carry tension. The first step in computing this value is to 
determine the distance from the bottom fiber to the neutral axis. This can be 
done using a ratio of the table stresses and the rib depth: 



yb = 



413 psi 

(413 + 1096) psi 



(32 in) = 8.75 in 



From this, the total tensile force is: 

T = 1 (413psi)(18in)(8.75 in) = 32,5001b 

This agrees exactly with the table value. With this value for tension, and 
assuming a yield stress in the steel of 50ksi, the required reinforcing steel area 
would be less than one square inch. The ACI code for minimum reinforcing 

steel area would supersede, thus requiring a reinforcing steel area of:"'® 

pmin = 0.005 = 

bd 

Solving for using the rib cross section yields: 

As = 0.005 (576 in2) = 2.88 in2 

This explains the table reference to minimum required reinforcing steel and 
agrees with Tedesko’s recommended steel area. The tension force and 
required steel areas are calculated similarly for the other cases. For case three, 
the required tension force is zero since the entire cross section is in 
compression at the crown. 

Tedesko computes a safety factor against buckling as the last table value 
for each cross section. The computations are based on Dischinger's formula:"*® 
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V = 33.21 El 
* Ha2 

where: H = Horizontal Thrust 

E = Effective Modulus of elasticity 
a = One-half the Span Length 
I = Moment of Inertia at the Crown 

Using the data for the first cross section with Tedesko’s thrust value of 499,820 lb 
yields a buckling safety factor of: 

Vs = 8.8 

This value compares very well with the table value of 8.7. The safety factors for 
the other cross sections are computed using this formula and show how the they 
are increased when the shell is moved to the rib extremity. Even with the 
reduced rib in case four, the safety factor against buckling is larger than when 
the shell is positioned at mid-height. 

Using the same hangar model and data as Whitney does, Tedesko has 
reached the opposite conclusion. In his table, with the shell moved to the rib 
extremity, the rib size can be reduced. Which conclusion is correct? 
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Chapter Four 
Resolving the Conflict 



Our examination thus far reveals questions which must be addressed to 
resolve the shell positioning conflict. In addition to these, other items which 
affect the issue are mentioned in the Whitney and Tedesko articles, but are not 
incorporated into their tables. 

First of all, Tedesko states that additional analysis is required: 

“. . . the writers do not believe it justified to assume that an investigation 
of the crown of the arch alone can determine the most economic cross 
section. Not only do the maximum moments vary in sign and magnitude 
along the arch axis, but also the relative importance of the volume 
change moments varies. In the lower quarters of the arch the volume 
change moments are only a small percentage of the total design 

moments.”^® 

Although this is a very serious discrepancy, he does not make the full span 
analysis he claims is necessary. 

Secondly, both sides mention arch deformation effects. And although 
they both imply that the deflections are easily approximated, neither side 
presents any relative data. Whitney even stresses the importance of 
investigating the moments caused by deformations, especially at the crown 

section where “the greatest increase in stress due to deflection occurs.”^"' 

Tedesko gives some justification for not investigating the additional moments: 

“. . . deformation moments are of important influence only for arches of 
small buckling safety and for arches which do not follow the pressure line 
for dead load” 

He obviously does not consider the deflection moments to be significant in this 
design. 
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Also, the two factions address using a reduced Young’s Modulus to 
account for concrete creep when computing deflections. Whitney states that the 
Young’s Modulus value for concrete should be reduced by two-thirds to three- 

quarters to calculate deflections under permanent load.^^ Tedesko suggests 
using a value of 2,000,000 psi for E^, to account for creep.^"^ 

To resolve the question of shell position, we will use an approach that 
incorporates analysis methods not available to Whitney and Tedesko and 
addresses the additional points mentioned above. Using the four step process 
outlined below, we will create a table for each of our three cross sections which 
we can use to make comparisons between the two design methods. The 
resulting tables are attached as appendix B. 

REVISED FOUR STEP METHOD OF ANALYSIS 



Step 1 Computina Forces. Moments and Deflections 

We make a full span analysis of the barrel shell roof section utilizing the 
Finite Element computer Program P-FRAME. From the finite element analysis, 
we are able to determine dead load moments, positive and negative moments 
due to different live load distributions, volume change moments due to 
temperature change, axial thrusts at each section and deflections due to the 
loadings. We model three different cross sections and create a table for each. 

We model the arches using the 20 foot width and 4 inch shell thickness 
Whitney specified. The first cross section has the shell at the mid-height of the 
rib and will be referred to as the “Whitney Arch”. The second is 14 feet long, has 
the shell located at the lower rib extremity and will be referred to as the 
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Tedesko Arch”. Both of these arches have ribs that vary from 1 8 / 32 at the 
crown to 18 x 40 at the springing. The third section is 14 foot v/ide, has a 
smaller rib and has the shell located at the lower extremity. It will be referred to 
as the “Reduced Tedesko Arch”, We chose to position the shell at the bottom of 
the rib for the Reduced Tedesko Arch since this will result in a compressive 
stress distribution through more of the arch span. Since v/e are designing in 
concrete, this is an important consideration. The rib for this section varies from 
9 X 32 at the crovm to 9 x 40 at the springing. The three cross sections are 
shown in Figure 7. 

We used 29 nodes to model each arch, 23 of v/hich are spaced 
horizontally from zero to 220 feet at equaJ 10 foot intervals. The other six nodes 
are placed to allow for the live loads necessary to make a full span analysis. To 
ensure symmetry, v/e placed three nodes on each side of the mid-span. Vertical 
positions for the nodes v/ere computed using the equation for a parabola. A 
one-line diagram of the model is shov/n in Figure 5. 

The moments of inertia at the crov/n for the first tv;o arches are taken 
directly from Whitney’s table. We use the corrected moment of inertia from the 
Tedesko table for the third. Moment of inertia at the springing is computed using 
Whitney’s dimensions. These calculations appear in Appendix A. To represent 
the varying cross section, v;e use a linear interpolation betv^een the crov/n and 
springing, adjusting the value every 10 feet. The Young’s f/odulus for all cases 

is 4 X 10^ psi, the same as used by Whitney and Tedesko, and the coefficient of 

thermal expansion v;e chose is 5.5 x 10'^ in/in/'F, the value that V/hitney 
recommends. The end restraints for the arches are modeled as fixed against 
both rotation and translation. For all three cases, v/e chose a linear elastic 
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analysis. 

To compute dead loads, we input the normal density of concrete, 150 
and P-FRAME computes the weights based on cross sectional areas and 
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Figure 5. Crcrwn cross sections and Ooe-Hne diagram for trve corrpoter rroOel. 
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lengths between the nodes. For the Tedesko arches, an additional externally 
applied dead weight was added to account for the reduced effective width. This 
was modeled as a uniformly distributed horizontal load. The values for this load 
were computed as: 



The live loads are modeled as uniformly distributed and externally 
applied using Whitney’s magnitude of 30 psf. For the full span analysis, we had 
to load the arch with several different live load distributions to produce 
maximum positive and negative moments at the points we wanted to 
investigate. Although data is available through P-FRAME for every nodal point, 
we concentrated our live load analysis on three significant points; the springing, 
the quarter point, and the crown. Each arch is loaded using known distributions 

to produce maximum positive and negative moments at the points of interest.^^ 
Load distributions used are shown in Figure 6. 

To model the volume change moments, we applied a uniform 
temperature change of -40 °F along the full span of each arch. The rib 
shortening contribution to the volume change moments is computed directly by 
P-FRAME. The program determines the axial deformations due to the applied 
dead and live loads. 

Step 2 Computing Deflection Moments 

To compute the additional moments due to deflections,we applied the 
theories of Robert S. Rowe and Milo S. Ketchum to the output data from P- 



(20 ft- 14 ft) (4 in 




yL)(i5oil2-| = 300i^ 
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FRAME. Both of these procedures are based on series approximations of the 
deflections. To account for the affect of creep, we used a Young’s modulus of 
2,000,000 psi in the computer analysis. 
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Figure 6. Live load distributions used for Arch analyses 

Ketchum’s procedure, published in the American Society of Civil 
Engineers (ASCE) Transactions, provides a method for computing final 
deflections as a function of initial deflection, moment and axial force.^s The 
derivation is based on a beam which is loaded both axially and laterally as 
shown in Figure 7. In the figure, is the moment due to q, the lateral load, wj 
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is the deflection caused by the lateral load, is the deflection due to the axial 
load and is the total deflection. Relating the deflections: 

Wo = Wj + Wa 

Assuming the elastic deflection curves for both loads are similar to their bending 

moment diagrams, a relationship is established between the ratios of the 

deflections and moments at the midpoint of the beam: 

Wa ^ P Wq 
wi Ml 



Solving for Wa: 



Wa 



Wj Pw, 
Ml 



Since Wa = w„ - wj; 



Wo 



w; + Wj 



PWo t 

Ml) 



Regrouping: 



Wo 




Ml) 



Wj 



Solving for 



Wo 







P^\ 

Ml) 



Thus, the final deflection, Wq, can be computed if the initial deflection, the axial 
force and the moment are known. This formula is adapted for our use in 
computing deflection moments for the tables in appendix B as follows: 
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The deflection moment can be computed if the initial deflection, axial force and 
moment at a section are known. The initial data (P, wj, MJ is available from our 
computer analysis. 




Figure 7. Axially atxJ Laterally loaded beam used in the Ketchum derivation. 



Robert S. Rowe’s procedure is derrived using a beam loaded both axially 
and laterally as shown in Figure 8. Rowe’s method expanded on Ketchum’s 

work by applying it to arches as well as beams.^^ Rowe’s procedure is once 
again based on the idea of a series of moments. From Figure 8, the deflection, 
Ab, due to the arbitrary lateral load can be expressed as: 




where “n” is a bending moment diagram shape factor. 



The moment caused by this deflection and the given axial load is: 



^«til ~ Pi 
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Figure 8. Loaded Beam used in the derivatton of Robert S. Rowe's method 



This moment causes an additional deflection of: 




Which, in turn, causes an additional moment: 

M. = P. A«i = P. (na 

Assuming the elastic curve maintains the same shape so the n^’s are similar, the 
total moment equation becomes: 
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iM = M), + Pa A b + Pa A b Hj 






E I 



E I 



Since 



£^i<i 
Ell , this becomes: 



P L^l 

M = Mb + Pa A b 1 1 + Oa I £ I ] order terms which go to zero 



Pa O 



Multiplying through by 



( I ' ) where A = Oa -§-r 

'1 -A / El 



M = 



Mb(l - A) + PaAb(l - A^) 
1 - A 



WithA<l,A^ goes to zero, therefore: 



^ Mb(l - A) 4- PaAb 
1 - A 



Assuming the axial load and lateral load bending moment diagrams are similar, 

A _ „ MbL^ 

— n p j 

"a ~ "b, and since v/e knov/ , the moment equation becomes: 



M = 



Mb 1-na- 



PaL^ 

El 



•-¥i , 



1 -n. 



PaL^ 

El 



P 

1 - n 

^ El 
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M = 



Mb 

PaL2 

Therefore; ^ * 

The total bending moment, including the deformation effects, can be computed 
in terms of the bending moment at the section, the axial load at the section and 
a bending moment diagram shape factor. This is very similar to Ketchum’s final 
result. 

Rowe applies this to curved beams. He presents a chart that relates the 
displacement ratios in straight and curved beams to the h/L value in the curved 
beam. From his chart, we see that for arches with rise to span ratios of less 
than 0.15, the deflection in an arch is less than 2% different from that in a 
straight beam.^® 



Step 3 Computing Cross Sectional Stresses 

Stress distributions are computed using the total moments and axial 
thrusts at each section and the standard P/A and My/I stress formulas. Results 
from these calculations appear in each of the tables. 

Step 4 Computing Buckling Safety Factors 

A revised method of computing buckling safety factors, based on the 
theories of Rowe and Ketchum, is presented in the following chapter. Data for 
each arch is presented at the end of this chapter. 
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EXPLANATION OF TABULAR DATA 

General 

We show a full span analysis in the three tables in appendix B. Data is 
computed for the crown, the left and right springing and the left and right quarter 
points. The data is displayed in columns from left to right along the arch length. 
An explanation of each line, along with a sample calculation for the “Whitney 
Arch” follows. The P-FRAME output file for the Whitney Arch is attached as 
Appendix C. 

Moment of Inertia. Line (A) 

Moments of inertia for the first two tables are computed based on 
Whitney’s data for the crown and springing rib sizes. The quarter point value is 
linearly interpolated. Data for the third table is calculated based on Tedesko’s 
reduced rib width, using Whitney’s variation in rib height . 

For Whitney’s Arch, the moment of inertia at the crown was previously 

computed as 58,000 in"^. For the given rib dimensions of 1 8 x 40 at the 
springing, the moment of inertia is 110,000 in^. Using a linear interpolation 
between the two to compute the quarter point value yields 84,000 in"^. 

Moments due to dead and live load. Line fBl 

These values are taken from P-FRAME output for the cross sections with 
positive moments acting clockwise at the left and counter-clockwise at the right 
hand end of a segment. The dead load is calculated based on a linear 
interpolation of cross section variation every 10 feet. Dead load for the full 20 
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foot width is applied for all three arches. Load positions for moments due to live 
loads are shown for the four cases in Figure 8. P-FRAME data is converted from 
ft-kips to in-lb for easy comparison with Whitney and Tedesko data and is 
rounded to four places. For the Whitney Arch at the crown (node 1 5): 



Moment due to dead load = -30.20 ft-kips = -363,000 in-lb 

Moment due to live load (1), 

(Maximum positive crown Moment) = -t-136.4 ft-kips = +1,636,000 in-lb 
Moment due to live load (2), 

(Maximum negative crown Moment) = -133.7 ft-kips = -1,604,000 in-lb 

Moment due to live load (3), 

(Maximum positive moment at the 
Left Quarter Point or Maximum 
negative moment at the Left 

Springing) = -64.48 ft-kips = -774,000 in-lb 

Moment due to live load (4) 

(Maximum negative moment at the 
Left Quarter Point or Maximum 
positive moment at the Left 

Springing) = +67. 10 ft-kips = +805,000 in-lb 



Initial Displacements. Line fCl 

These values are also taken from P-FRAME output and represent 
displacement of the nodal points from their initial positions due to the indicated 
loads. The displacements shown are in inches and do not include deformation 
moment effects. The dead load displacements are computed using E = 
2,000,000 psi to account for creep. For Whitney’s Arch at the crown: 



Initial Displacement due to Dead Load 
Initial Displacement due to Live Load case (1) 
Initial Displacement due to Live Load case (2) 
Initial Displacement due to Live Load case (3) 



= -0.345 inches 
= -0.541 inches 
= -1-0.419 inches 
= -t-0.180 inches 
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Initial Displacement due to Live Load case (4) = -0.302 inches 



Axial Thrusts. Line (D) 

Once again, the Values are taken directly from P-FRAME output. The 

thrusts are normal to the indicated cross section and are listed in kips. For 

Whitney’s Arch at the crown: 

Axial thrust for Dead Load =351.6 kips 

Axial thrust for Live load Case (1 ) = 60.6 kips 

Axial thrust for Live load Case (2) = 70.6 kips 

Axial thrust for Live load Case (3) = 41.1 kips 

Axial thrust for Live load Case (4) = 90.1 kips 

This would give a total thrust of 412,200 lb for dead plus live load for maximum 
positive moment at the crown, and a thrust of 482,800 lb for dead plus full span 
live load. These values are 3% higher than the thrusts of 400,400 lb and 
469,600 lb computed using Whitney's formulas. 



Temperature Change Moments. Line (El 

These moments are computed based on the 40°F temperature drop used 
by Whitney applied over the entire arch. The rib shortening contribution, 
computed by Whitney for his volume change moments, is not included here 
since it is computed as part of the Load Moment in line (B). The P-FRAME 
output is once again adjusted from ft-kips to in-lb for comparison with Whitney 
and Tedesko table values. For the crown section of Whitney’s Arch; 

Temperature Change Moment = 59.87 ft-kips = 718,000 in-lb 

If our computed values of dead load, live load for maximum moment at 
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the crown, and temperature change moments are added together, we can make 
a comparison with the total moment value listed in Whitney’s table. The P- 
FRAME total moment would be 1,991,000 in-lb. This is 27% lower than 
Whitney’s table value of 2,736,000 in-lb. 

Moments due to Deflection. Line (F) 

As the h/L value for our arch is 0.125, according to Rowe, we can 
compute the deflection moments using Ketchum’s method. We use P-FRAME 
output to compute the axial thrust, initial deflection and moment at the section for 
each live load condition. The dead and temperature change loads will not 
create deflection moments since they are uniformly distributed across the entire 
arch span. The arch axis will not deflect from the funnicular line under these two 
loads. 

For Whitney’s Arch at the crown, the deflection moment at the crown due 
to live load case (1) is computed using; 



From P-RFRAME output 

for Live load case (1) P= 60,600 lbs (compression) 

Wj = -0.541 in 
Ml= 1,636,000 in-lb 



The deflection moment at the crown for live load case (1) is: 

M = + 33,000 in-lbs 

The moment sign is determined by the orientation of the deflection and axial 
force. In this case, the arch is deflecting downward, and the axial force is 
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compressive, thus, a positive moment results. 



Worst Case Moments (G) 

The worst case moments are totaled from Lines B, E and F for dead and 
live loads. The live load case which produces the largest appropriate moment 
is used in each computation. A circled number appears next to the worst case 
moment value to indicate which live load was used. For the crown section of 
Whitney’s Arch; 

Worst case total Moment = +2,024,000 in-lb 

The worst case total moment for the crown results from a live load for maximum 
positive moment. 

Section Stress. Line (H) 

Stresses are calculated based on worst case moments and associated 
axial loads. Standard stress formulas are used with the computed cross 
sectional areas and section moduli for the appropriate point in the arch. For the 
crown section of Whitney’s Arch: 

Axial Load for Worst case Total Moment = 412,200 lb (compression) 

Top and bottom fiber stresses: 

fj = - 840 psi and fb ~ 

Buckling Safety Factor 

A safety factor against buckling is computed based on the Initial 
Deflection Method for arches developed in Chapter 5. First, models of the 
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Whitney Arch, the Tedesko Arch and the Reduced Tedesko Arch are given a 
parabolic imperfection. We determine the critical load using the Initial 
Deflection Method, and compute a buckling safety factor by comparing the 
critical load to the dead load plus the 600 Ib/ft live load used by Whitney. The 
finite element analysis program P-FRAME is used to compute the initial 
deflections. 

For each arch, the initial imperfection is a parabola with a midspan rise of 
0.625 ft. We compute the revised nodal coordinates using the transformation 
described in Chapter 5. The revised coordinates are input to the P-FRAME 
program, and we apply a uniformly distributed horizontal load across the entire 
span. To compute the critical load, we must eliminate the rib shortening 
contribution. A purely parabolic model of the arch is loaded with the same 
uniform load, and the resulting deflections are subtracted from those computed 
using the offset model. The final deflections are due to bending moment only. 
The critical load occurs when the deflection is equal to the offset. 

For Whitney’s Arch, at a uniform load of 23.5 kips/ft, the deflections for the 
offset model at the three-quarter point (node 23) are: 

6 j (23 = -0.224 in and 5y23 = -2.49 in 

The deflections for the parabolic model are: 

5x23 = -2.76 in and 5y23 = -9.56 in 

This results in a total deflection of: 

5 = -7.51 in =-0.625 ft 

To compute the buckling safety factor, we assume a uniformly distributed dead 
load equal to the average cross sectional value. The buckling safety factor for 
Whitney’s arch is: 
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Safety Factor = 23.5 /2.2 = 10.7 

The P-FRAME input and putput used in this calculation are attached as 
appendix D. Using the Initial Deflection Method for the other two arches yields 
critical loads of: 

Tedesko Arch = 47.5 kips/ft 

Reduced Tedesko Arch = 27.2 kips/ft 

The resulting Buckling Safety Factors are: 

Tedesko Arch = 25.0 

Reduced Tedesko Arch = 17.4 

These values are higher than the Dischinger values listed in Tedesko’s table, 
but show the same trend. Even with a reduced rib, the Tedesko cross section 
has a greater safety against buckling than Whitney’s. 



37 



Chapter Five 



The Initial Deflection Method 
For Computing Critical Buckling Loads 



Using the ideas of Rowe and Ketchum, we will develop a relationship 
between successive deflections and buckling. The general formula is 
developed using an axially loaded column. 

An axially loaded column which is perfectly straight theoretically will not 
buckle, regardless of the applied load. It will only deform along its axis in 
accordance with the well known formula: 



where: P = Axially applied load 

L= Column length 
A = Cross sectional area 
E = Young’s modulus 

Given some type of initial imperfection, however, the column will buckle under 
sufficient load. For the axially loaded column shown in Figure 11 that has an 

initial parabolic offset with a value of 5 q at mid-height, the initial moment at the 
mid-span is: 



Mo = P5o 

This moment , in turn, causes an additional deflection, 6^: 

5 = MqL^ ^ (P 5q) 

^ 12EI 12EI 
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Letting 



B = 

12 E I, the additional moment due to the deflection 5-| is: 



Mi= P5i = P5oB 



? 




? 

Figure 9. Axially loaded column with initial parabolic offset 



This moment will create an additional deflection , Sj: 



52 = 



Ml 

12 El 



6oB2 



The total moment is: 

Mt = M 0 +M 1 + M 2 + . . . = P5o+P5i + P52 + . . . = P5o (1 + B + B^ + . . . ) 

Thus, if B is greater than one, the series diverges, the moments will grow without 
bound, and the column will buckle. The critical load occurs when B = 1 . When 

B = 1, the original offset, 5o, and the initial deflection, 5j, will be related as: 
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5i = B 5o = 5o 



We have defined the critical point in terms of deflection. When the initial 
deflection is equal to the original offset, the column is at the critical point. If the 
initial deflection is greater than the original offset, the column will buckle. This 
criteria is easily applied to output from finite element computer programs. To 
check buckling, one only needs to compare the computer generated deflections 
to chosen input imperfections. 

As a check, we will compare the results from the “Initial Deflection 
Method” with several well known buckling formulations. 



Euler Column Buckling 

We will check the initial deflection buckling criteria against the classic buckling 
problem presented by Euler. His solution for the critical load of a column 

hinged at both ends is; 

p _ jr^EJL 

'' L2 

where: P = Column axial load 

L = Length between the supports 
E = Young’s Modulus 

I = Minimum moment of inertia of the cross section 
The physical model used is a 100 inch long steel beam with a 12 square inch 
cross section which is 6” x 2 “. The beam is hinged at both ends. 

Since the chosen cross section results in a minimum moment of inertia of 

4 in^, the critical load using Euler’s formula is: 
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Per = 118 kips 



To test the Initial Deflection Method, the beam is modeled on the finite 
element computer program SAP-90 using frame elements with 11 equally 
spaced vertical nodes. The horizontal offset for each node is calculated using a 
parabolic equation with the chosen maximum offset as 0.625 inches at the mid- 
span. The column is modeled with a pin at the bottom and a roller at the top. 
Loads are applied in the negative vertical direction at the top node. The critical 
load is one that produces a 0.625 inch deflection at the middle node. To 
determine the critical point, we start at the theoretical critical load and perform 
iterations until we read a deflection of 0.625 inches on the output. 

Using the Initial Deflection Method with SAP-90, the critical load is : 



Next, we check the initial deflection buckling criteria against the theory for a 
simply supported plate. The physical structure we model is a 100 inch by 100 
inch steel plate which is one inch thick. 

The general formula for critical load per length for a simply supported 
plate uniformly compressed in one direction is:^*^ 



Per = 116 kips 



This is 1 .7% lower than Euler’s theoretical value. 



Plate Buckling 




where: 



a= Horizontal plate dimension 
b = Vertical plate dimension 

= Load per unit length along the horizontal 
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and; 



where: E = Young’s modulus 

h = Plate thickness 

V = Poisson’s ratio 

For a square plate, a = b, and this becomes: 

(Nj„ = 

For the chosen plate, with E = 30 x 10^ psi and v = 0.3, the flexural rigidity is: 

D = 2,747 in-kips 

From the Timoshenko formula, the critical distributed load is: 

(N,)c= 1084 

The structure is modeled using the finite element computer program SAP-90 
with 100 plate elements. The computer model is composed of 121 nodes with 
11 nodes spaced equally along the horizontal, and 11 vertical nodes equally 
spaced at each of these. Each vertical line of nodes has a parabolic offset with 
a maximum of 0.625 inches at the middle node. The plate is simply supported 
on all sides with a pin along the bottom edge and rollers along the other three. 
We model the uniform load across the top of the plate as a pressure load. 

As with the axially loaded column, the critical loading occurs when the 
mid-span horizontal deflection equals the original offset. Using this criteria, we 
load the plate at the theoretical critical load, and adjust as necessary until we 
see a mid-span deflection of 0.625 inches in the computer output. 

Using the Initial Deflection Method, the critical load is: 
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(N,)ct= 911 ^ 

in 

This result is 15.6 % lower than the theoretical buckling load. This result is 
suprising in light of the success with the column data. 
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Arch Buckling 



In a two hinged parabolic arch, uniformly distributed loads are carried 
axially to the supports. Bending moments are zero throughout, and the arch 
simply “squats” due to rib shortening. By introducing an initial imperfection as 
we did with the axially loaded column, however, bending moments are created 
(Figure 1 2 refers). The moments can be expressed as: 

Mo= P5o 

where: P = Axial load at the section 

5q = initial offset from the funicular line 

The greatest moments will occur at the points where the largest offset is located. 
These bending moments will, in turn, cause additional deflections. These 
deflections are expressed as they were for the axially loaded column as: 

8 = MqL^ = (p5q)l^ 

^ 12EI 12EI 

Using the same derivation as in the case of the axially loaded column, we can 
define the critical point of arch buckling in terms of deflection. When the initial 
deflection is equal to the original offset, the arch is critically loaded. If the initial 
deflection is greater than the original offset, the arch will buckle. To compute 
buckling loads using the Initial Deflection Method, we need to compare 
computer generated deflections to our chosen offsets. 

To validate the Initial Deflection Method for arch buckling, we will 
compare it with Timoshenko’s theoretical results. The model we use is a 2 
hinged steel arch with a constant 6” x 2” cross section. The arch spans 100 
inches horizontally and has a rise of 10 inches. 
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Timoshenko’s general formula to compute critical loads for a uniformly 



loaded parabolic arch with a constant cross section is:^’ 

qcT = ^ ^ 



where: E = Young’s modulus 

I = Moment of inertia 
L= Horizontal span length 

>.4 = A factor depending on the height to span 
ratio and the number of hinges 
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Figure 10, Bending Moments in Arches caused by Inrtial Offset 
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From Timoshenko’s Table 7*5 with h/L = 0.1, = 28.5. The critical load with E 

= 30 X 10^ is; 

q cr = 3.42 kips per inch of horizontal span 

The Initial Deflection Method is tested using the finite element computer 
program P-FRAME. We model the arch using 21 nodes equally spaced along 
the horizontal and compute initial vertical coordinates for these nodes using a 
parabolic equation with a rise of 10 inches at the middle node. The arch is free 
to rotate and restricted from translating at both ends. 

We will apply a parabolic offset to each half span as shown in Figure 11. 
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Figure 11. Parabolic Offset used in the Computer Model 

The maximum offset will be 0.625 inches at each quarter point. Since we want 
the nodes to be offset from the initial parabolic curve, each node will have to be 
adjusted both vertically and horizontally as shown for node 5 in Figure 12. The 
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change in coordinates will be a function of the parabolic offset, 85 , and the angle 



65 as: 



5yS = 5s (cos 65 ) 



and 



^5-^5 (jj,^ 0^1 




Figure 1 2. Blown up view of the coordinate transformation at node 5 

The arch is loaded using uniformly distributed horizontal loads. To 
eliminate the affect of rib shortening from the computer output, a purely 
parabolic arch model without initial offsets is loaded with a uniformly distributed 
critical load. The rib shortening deflections are subtracted from the deflections 
computed using the offset model. The resulting deflections are due to bending 
moment only. Once again, the theoretical buckling load is chosen as a starting 
point and iterations are performed until we see an adjusted deflection of 0.625 
in either half of the arch. Using the Initial Deflection Method, the buckling load 
for the arch model is: 
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Pe, = 2.92 ^ 
in 

This is 14.6% lower than Timoshenko’s theoretical value of 3.42 kips/in. 
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Chapter Six 
Conclusions 



From the debate between Charles S. Whitney and Anton Tedesko we 
can draw several conclusions. First of all, Whitney’s 1925 article was a good 
guideline for arch design. His formulas and charts give information which is 
backed up by modern methods of analysis. Using only graphic statics and the 
Calculus, he computed formulas and charts to determine thrusts and moments 
for significant loads and load positions at several points in the arch. For 
maximum positive moment at the crown, his distribution is less than 1% different 
from accepted design guidelines used from the 1950’s to today. His adjustment 
factors allow quick computations for a very diverse range of arch designs. 
Although his conclusions for thin shell design are vehemently challenged by 
Tedesko, his method of calculating arch thrusts, live load moments and volume 
change moments are not. 

Secondly, Tedesko’s claim that a full span analysis is necessary for arch 
design is valid. By looking at our full span analysis, we can see that stress 
distributions at the springing, quarter point and crown must all be investigated. 
The crown analysis that Whitney and Tedesko both present is not sufficient to 
design a barrel shell roof. 

Next, both men were correct in their claim that the deflection moments 
are not significant for Whitney’s hangar model at the assumed live loads. We 
can see, however, that the deflection moments and corresponding stresses 
must be addressed in barrel shell roof design. As the factor of safety against 
buckling is reduced, the deflection moments become more significant. The 
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methods of Rowe and Ketchum are a reliable way to determine deflection 
moments from data available from most finite element programs. 

Additionally, Whitney’s emphasis on volume change moments is valid. 

His claim that the shell should be located at the mid-height of the rib to reduce 
volume change moments is also valid. There is an irrefuteable relationship 
between the moment of inertia and the volume change moment. When we look 
at stress distributions over the full span and take all loads into account, 
however, we see that the decreased moment of inertia has a drawback. 
Reducing the moment of inertia increases the tensile stresses near the 
supports. The tensile stress at the springing for Whitney’s Arch is double that of 
Tedesko’s Arch. 

Also, there is a significant advantage in placing the shell at the bottom of 
the rib as opposed to the top. Because concrete is weak in tension and strong 
in compression, we want to reduce tensile stresses as much as possible. 
Placing the shell at the bottom of the rib does this. We can see from the stress 
distributions for Tedesko’s Arch and Tedesko’s Reduced Arch that tension only 
exists in upper part of the cross section for the first and last quarter of the arch. 

The Initial Deflection Method is a tool which can be used with existing 
finite element computer programs to compute buckling safety factors for a 
variety of structures. Although longhand computations of initial offset values are 
tedious, this could easily be written into a finite element computer program. The 
critical loads computed with the Initial Deflection Method are conservative for 
plates, but could provide a ready check for buckling capacity for arches during 
an iterative design process. 
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APPENDICIES 



APPENDIX A: CALCULATIONS FOR MOMENTS OF INERTIA 

APPENDIX B: DATA TABLES FOR THE FOUR STEP METHOD OF ANALYSIS 

APPENDIX C: SAMPLE COMPUTER DATA FILE FOR WHITNEY’S ARCH 

APPENDIX D: SAMPLE COMPUTER DATA FILE FOR INITIAL DEFLECTION 

METHOD CALCULATIONS 



Moment of Inertia Calculations 



Section 1 : 18" thick rib with shell located at mid-height 



As = As’ = 3.0 inches, n = 7 





lo 


A 


7t 


Ayt 


d to na 


Ad^ 


Rib 


49,200 


576 


16 


9220 


0 


0 


Left SheU 


592 


444 


16 


7100 


0 


0 


Right Shell 


592 


444 


16 


7100 


0 


0 


Top Steel 


2 


21 


2.5 


53 


13.5 


3830 


Bottom Steel 


2 


21 


29.5 


620 


13.5 


3830 






y = 16.0 in 




1 = 58,000 in^ 







Section 2: 18" thick rib with shell located at top (Range 70% effective) 



As = As’ = 


3.0 inches. 


n = 7 










lo 


A 


7t 


Ay^ d to na 


Ad2 


Rib 

and Steel 


56,800 


618 


16 


9890 6.9 


29,400 


Left SheU 


400 


300 


2 


600 7.1 


15,100 


Right SheU 


400 


300 


2 


600 7.1 


15,100 






y = 9.1 in 




1= 117,200 in^ 
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Section 3: 27" thick rib with shell located at too (Flange 70% effective) 



As = As’ = 4.5 inches, n = 7 





lo 


A 


7t 


Ayt 


d to na 


Ad^ 


Rib 


73,700 


864 


16 


13,800 


5.3 


24,300 


Left Shell 


376 


282 


2 


564 


8.7 


21,300 


Right Shell 


376 


282 


2 


564 


8.7 


21,300 


Top Steel 


3 


32 


2.5 


80 


8.2 


2,150 


Bottom Steel 


3 


32 


29.5 


944 


18.8 


11,300 






y = 10.7 in 




I = 154,800 in’^ 







Section 2: 9" thick rib with shell located at top (Flange 70% effective) 



As = As’ = 


1.5 inches, n 


= 7 












lo 


A 


7t 


Ayt 


d to na 


Ad2 


Rib 


24,580 


288 


16 


4608 


9.6 


26,740 


Left SheU 


450 


318 


30 


9540 


4.4 


6160 


Right Shell 


450 


318 


30 


9540 


4.4 


6160 


Top Steel 


0 


1.5 


2 


3 


21.6 


700 


Bottom Steel 0 


1.5 


30 


15 


4.4 


30 



y = 25.6 in I = 65,300 
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Whitney Arch 





Uft 

Springing 


Uft 
Qtr Pt 


Crown 


Right Right 

Qtr Pt Springing 


A) Moment of Inertia 


110,000 


84,000 


58,000 


84,000 


110,000 


(in4) 

B) Load Moment (xio^ in-ib) 


Dead Load 


-1.936 


-bO.321 


-0.363 


+0.321 


-1.936 


Live Load 1 


-b2.575 


-1 .535 


+1.636 


-1.535 


+2.576 


Live Load 2 


-2.812 


+1 .492 


-1 .604 


+1.492 


-2.812 


Live Load 3 


-6.572 


+2.996 


-0.774 


-1.686 


+4.159 


Live Load 4 


-b6.335 


-3.010 


+0.805 


+1.643 


-4.396 


C) Initial Displacements (in) 


Dead Load 


0 


-0.476 


-0.345 


-0.476 


0 


Live Load 1 


0 


+0.196 


-0.541 


+0.196 


0 


Live Load 2 


0 


-0.260 


+0.491 


-0.259 


0 


Live Load 3 


0 


+1.155 


+0.180 


-0.917 


0 


Live Load 4 


0 


+1.092 


-0.302 


+0.980 


0 


D) Axial Thrust (kips) 


Dead Load 


397.4 


362.3 


351.6 


362.3 


397.4 


Live Load 1 


61.8 


62.8 


60.6 


62.8 


61.8 


Live Load 2 


85.0 


72.5 


70.6 


72.5 


85.0 


Live Load 3 


56.8 


42.9 


41.1 


41.2 


39.7 


Live Load 4 


89.9 


92.4 


90.1 


94.0 


107.1 


E) Temperature Change 


-1 .774 


+0.096 


+0.718 


+0.096 


-1 .774 


Moment (xio^ in-ib) 
F) Deflection Moments 


Live Load 1 


0 


-0.012 


+0.033 


-0.012 


0 


Live Load 2 


0 


+0.019 


-0.030 


+0.019 


0 


Live Load 3 


0 


+0.050 


-0.007 


-0.039 


0 


Live Load 4 


0 


-0.104 


+0.028 


+0.098 


0 


G) Worst Case Total 


Moment (xio^ in-lb) 


10.38 (3) 


+3.433(3) 


+2.024(1) +2.158(4) -8.106(4) 


H)Section Stress (psi) 


Top 


-b1593 


-996 


-840 -755 


+1168 


Bottom 


-2144 


+475 


+277 +169 


-1780 
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Tedesko Arch 





Uft 

Springing 


Left 

QtrPt 


Crown 


Right 

QtrPt 


Right 

Springing 


A) Moment of Inertia 


218,400 


167,400 


116,300 


167,400 


218,400 


(in4) 

B) Load Moment (xio^ in-ib) 


Dead Load 


-2.354 


+0.268 


+0.002 


+0.268 


-2.354 


Live Load 1 


+2.459 


-1.527 


+ 1.687 


-1.527 


-2.459 


Live Load 2 


-2.941 


+1.497 


-1.555 


+1.497 


+2.941 


Live Load 3 


-6.643 


+2.972 


-0.744 


-1.685 


+4.080 


Live Load 4 


+6.160 


-3.002 


+0.876 


+1.655 


-4.562 


C) Initial Displacements (in) 


Dead Load 


0 


-0.474 


-0.663 


-0.474 


0 


Live Load 1 


0 


+0.078 


-0.309 


+0.078 


0 


Live Load 2 


0 


-0.155 


+0.165 


-0.155 


0 


Live Load 3 


0 


-0.595 


+0.064 


+0.446 


0 


Live Load 4 


0 


+0.518 


-0.208 


-0.524 


0 


D) Axial Thrust (kips) 


Dead Load 


400.5 


365.5 


354.5 


365.5 


400.5 


Live Load 1 


61.3 


62.3 


60.1 


62.3 


61.3 


Live Load 2 


84.5 


71.9 


70.0 


71.9 


84.5 


Live Load 3 


56.6 


42.6 


40.7 


40.9 


39.4 


Live Load 4 


89.3 


91.7 


89.4 


93.3 


106.4 


E) Temperature Change 


-3.506 


+0.193 


+1.424 


+0.193 


-3.506 


Moment (xio^ in-lb) 
F) Deflection Moments 


Live Load 1 


0 


-0.005 


+0.019 


-0.005 


0 


Live Load 2 


0 


+0.011 


-0.012 


+0.011 


0 


Live Load 3 


0 


+0.026 


-0.003 


-0.018 


0 


Live Load 4 


0 


-0.048 


+0.019 


+0.050 


0 


G) Worst Case Total 


Moment (xio^ in-lb) 


-12.50 (3) 


+3.481(3) 


+3.132(1) +2.166(4) -9.939(4) 


H)Section Stress (psi) 


Top 


+1261 


-844 


-958 


-684 


+897 


Bottom 


-1028 


-96 


-95 


-219 


-923 
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Reduced Tedesko Arch 





Left 

Springing 


Left 
Qtr Pt 


Crown 


Right 
Qtr Pt 


Right 

Springing 


A) Moment of Inertia 


123,800 


94,600 


65,300 


94,600 


123,800 


(in4) 

B) Lx>ad Moment (xio^ in-ib) 


Dead Load 


-1.506 


+0.135 


+0.001 


+0.135 


-1.505 


Live Load 1 


+2.509 


-1.532 


+1.663 


-1.532 


-2.509 


Live Load 2 


-2.891 


+1.496 


-1.573 


+1.496 


+2.891 


Live Load 3 


-6.617 


+2.970 


-0.755 


-1.684 


+4.113 


Live Load 4 


+6.235 


-3.005 


+0.846 


+1 .649 


-4.495 


C) Initial Displacements (in) 


Dead Load 


0 


-0.507 


-0.715 


-0.507 


0 


Live Load 1 


0 


+0.153 


-0.520 


+0.153 


0 


Live Load 2 


0 


-0.256 


+0.325 


-0.256 


0 


Live Load 3 


0 


-1.041 


+0.132 


+0.800 


0 


Live Load 4 


0 


+0.938 


-0.327 


-0.903 


0 


D) Axial Thrust (kips) 


Dead Load 


318.9 


291.7 


282.9 


291.7 


318.9 


Live Load 1 


61.5 


62.5 


60.3 


62.5 


61.5 


Live Load 2 


84.7 


72.2 


70.2 


72.1 


84.7 


Live Load 3 


56.7 


42.7 


40.9 


41.0 


39.6 


Live Load 4 


89.5 


92.0 


89.7 


93.6 


106.7 


E) Temperature Change 


-1.989 


+0.107 


+0.805 


+0.107 


-1.989 


Moment (xio^ in-ib) 
F) Deflection Moments 


Live Load 1 


0 


-0.010 


+0.032 


-0.010 


0 


Live Load 2 


0 


+0.019 


-0.023 


+0.019 


0 


Live Load 3 


0 


+0.045 


-0.005 


-0.033 


0 


Live Load 4 


0 


-0.089 


+0.030 


+0.089 


0 


G) Worst Case Total 


Moment (xio^ in-lb) 


-10.11 (3) 


+3.257(3) 


+2.501 


(1) +1.980(4) -7.989(4) 



H)Secrion Stress (psi) 

Top +2189 -1355 -1350 -1001 +1600 

Bottom -1074 -94 -126 -245 -978 
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< ‘-T inter Ifanner bc'e= Here ■ 
LIVE LOAD ANALYSIS (4 LOAD CASES) 
l-iHITNEY'S AF'<Ch 



Case Resul ts 




Joint 


Load 


X-Di sp 1 . 


Number 


Case 


(in) 


1 1 


1 


- . 02 i BO 




i 


. 02769 




3 




IE 


1 


-.011 07 




E 


•01>17 






. 353d6 






-.3^SE5 


13 


1 


- . OOGO^i 




c. 


. 0 1 23 1 






. 3^6‘-ro 






-.3^169 


1^ 


1 


. 000 1 3 






. 00250 

H r' 



-.31637 



15 


1 


O , CiriC;*’)! ) 
0 . 00000 
.309^9 




p 




2 


- . .6 '1)9 6 


1 6 


1 


- . C)00 1 ?; 




i... 


- . OoEE-O 




-• 


.21631 




■-t- 


-.31B9-f 


^ •“} 


i 


. ■'"•OGOr. 




is 


--.01281 




d 


. 3 3.3 7 B 




‘-t 


- . 33855 


IE. 


J 


. 0 1 1 ■:)? 






“ . 0 1 6 1 7 




d 


- . w4tiC; Cl 






1 




■ og^e'2 




•Z. 


. ‘.'CL , 6 




3 


^ 39739 




*■ '■ 


' . w’»jd2' .' 


E ‘ 


c. 


r L‘279'6 
-- . 03376 




E 


35 M 0 




■^1- 


- . 3573 J 


21 


1 


.05503 



:PiME-' Lxnf^ai'r Elastic analysis results 

;lle:y 



Y-Displ . 
( in ) 


Rotat ion 
( rad ) 


-.21379 
. 1 1537 
-.39907 
. 79570 


- . 00 1 5i:> 
. C'O 1 L8 
. CH;)6-+b 
- . 00255 


- . 32009 
.21063 

. 56238 


- . 00 1 
. 00 1 38 
. 00289 
- . 00298 


-.35952 
.29506 
— . 600' 30 

. 1 


- . 00 190 
. 00 1 32 

. ■')03C, 1 

- . ■■.)03 i 0 


- . 99092 
. 'd'7 1 7-0 
-.21'i7B 


- .. 00083 
. 'jOOTS 

. U'.'wr.'r^o 


-.59159 
. 9 1 ^-^60 
. 17959 
-.30152 


- V YQOp 
< ) . )■ )( ( ) 

. '.:.'03 1 6 
- . ;'0316 


- . 99092 
.3717^- 
.52219 
- . 6h737 


. 00083 

- - 1)0073 

- “ (>0255 


-.35952 
. 29306 
.7^3-r6 
- . 90992 


. 00 1 9^) 

- . I'C- 1 32 
. 0'“) 1 S*‘) 
-.001 71 


- . 3E0'j9 
.21063 
.,83958 
.9^399 


, or* 197 
- .00138 
. 00 1 63 
-.00159 


— . 1. 1 S7"i- 
. 1 153^:^ 
.926 j 1 
- i. - 0299":> 


. 00153 
-. OO 1 96 
. 00 1 1 9 

- 00 1 C'5 


■1 -1 c:. 

— . .. 1 CL 

, .j7.)'7v 
. 95o6i.‘5 

-1 ,0578£ 


. 00158 

- . ^ )'.> 1 9 - 
. .."n':.'91 

— . Ou('G 1 


. 00995 


. 00 i 37 
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F r 1 n t e r E- a r*i n e r b c< e = He r e 
LIVE LOAD ANALYSIS (A LOAD CAbELE ^ 



;.JHITNLY‘5 AF^'CH 



Zase Resul ts 
Joint Lb"ad 

vJumber Case 



3 

L 



c: 

3 



33 



2h 



1 

E 

3 

A 



X-Di sp 1 . 
( in) 

“ : 35931 

— . 3 c 6 7 



, 0833A 
03865 
,3A971 
, 35602 

, 09293 
, 09866 

z 1 ^ 



. I GOES 
- . 10531 

.3169-h 

- . dS300 

. 1023S 
10570 
. £6077 
~ . £6‘+OQ 



“D i sp 1 . 
(in) 


Ro tat ion 
( rad ) 


- . 09668 
1 .01375 
1 - 10099 


“ ■ 00 1 25 
- 000 OS 


. 19889 
~ . E205E 
. <^'‘7 1 05 
1 .092(b9 


1 1 

OOOO 


. 196':)0 
- . 98035 


. 00065 
- . 0005c: 
- . 1 06 
. 00 1 i 9 


.SE698 
- . EoE'IE 
. 3 “LEY 
- . 89977 


. 00036 

- . Oi'iOFA 

- . 1 35 

.00198 


.E3E1 1 
- . £78 1 6 


- .00016 
. 000c!.F 



, 6566V 
I 6 V 6 7 '+ 



cc 


i 


. 08602 


. 1 9“fL*U 


-- . (H'.055 




z' 


- .08730 


- . 2 i ^'’52 
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6 
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■ 


- . i 89 




*7 


- .. 13697 


■ . -j-oAbS 


.oOj 95 


rl / 


1 


. 05956 
- . 05606 


. 11607 


- „ 0C'07 1 
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— . i £ V c*: 5 


- 00079 




3 


. 10316 


. 2279 1 


2881^6 






10267 


-.2-7 109 


C. ■- 


y 




.03715 


- u 'JOC‘55 




z 


■■ . U 1 7bf j 
. 03622 


- . 09 1 36 
. 9561 

— . ^ 6 V '■vJ 6 ’ 


. C’O JciO 
. ‘J ' i U c. 




-i 


-.03098 


. 00107 






0 . (^0(‘C3> 


0 - C' 00^>0 


0 . (jOO('K‘ 
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O . ^>ooou 


0 . 00000 
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0 r 0 ‘ 1 
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0 . C> 0 C )00 
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Printer Banner boee Here 
LIVE LOAD ANALYSIS (A LOAD CASES) 
LHIINEY’S ARCH 



Case Results 
"Jo 1 n t Load 

Number Case 



1 



1 

id 

3 



1 

a 

d 



SUPPORT REACTIONS 



X-Reac t ion 
(kips) 

60.59S 
70 . 59 1 
90.93 
90 - 259 

-60 . 59B 
-70.591 
— 1 ^ 0 ^ 9 3 0 

r , {■ • cr 



Y-Reac t i on 
(kips) 

16. 5<Jt) 
99 . 500 

9b. 305 
1 9 . 695 



1 ^ 



99 





Z-Reac t ion 
(K-f t ) 

-E19.66E 
d39 . j63 
597.655 
-527.955 

d 1 9 .663 
—339 . 663 
396 .615 
—666.316 



^A^1E Linea'“ Elastic analysis results 
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Str 

05 Sep ?1 



No . 0 
1 :E8 p 



03 




I 



CRITICAL LOAD ON WHITNEY ARCH 



♦♦♦ MEMBER LOAD DATA *** 



’■ase 1 - 


member 


distributed 


loads 










Sloped UDL 


Pro j . UDL 


Local UDL 


Local UDL 


Triangular 


Thermal 


1 « 


K/ft 


slope 


K/ft horiz 


k/ft perp 


K/ft parll 


K/ft @ GJ 


Change ( F ) 


1 


0 




-23.5 


0 


0 


0 


0 


2 


0 




-23 . 5 


0 


0 


0 


0 


3 


0 




-23.5 


0 


0 


0 


0 


4 


0 




-23.5 


0 


0 


0 


0 


5 


0 




-23.5 


0 


0 


0 


0 


6 


0 




-23.5 


0 


0 


0 


0 


7 


0 




-23.5 


0 


0 


0 


0 


8 


0 




-23.5 


0 


0 


0 


0 


9 


0 




-23.5 


0 


0 


0 


0 


.0 


0 




-23.5 


0 


0 


0 


0 


.1 


0 




-23.5 


0 


0 


0 


0 


.2 


0 




-23.5 


0 


0 


0 


0 


.3 


0 




-23.5 


0 


0 


0 


0 


.4 


0 




-23.5 


0 


0 


0 


0 


.5 


0 




-23.5 


0 


0 


0 


0 


,6 


0 




-23.5 


0 


0 


0 


0 : 


.7 


0 




-23.5 


0 


0 


0 


0 


.8 


0 




-23 . 5 


0 


0 


0 


0 


.9 


0 




-23.5 


0 


0 


0 


0 


:o 


0 




-23.5 


0 


0 


0 


0 


:i 


0 




-23.5 


0 


0 


0 


0 


12 


0 




-23.5 


0 


0 


0 


0 


13 


0 




-23.5 


0 


0 


0 


0 


:4 


0 




-23.5 


0 


0 


0 


0 


:5 


0 




-23.5 


0 


0 


0 


0 


:6 


0 




-23.5 


0 


0 


0 


0 


17 


0 




-23.5 


0 


0 


0 


0 


18 


0 




-23.5 


0 


0 


0 


0 



>ped UDL, Projected UDL & Point Loads act in the global coordinate system. 
:al Perpendicular, Local Parallel, Triangular Loads act in 
i local member coordinate system. 

.angular Loads are 0 at the lower joint with the magnitude specified at 
; greater joint. 
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CRITICAL LOAD ON PARABOLIC ARCH 
WHITNEY’S ARCH 



JOINT DISPLAC: 






?s ♦♦♦ 



ase Results 
Joint Load 
iJumber Case 



X-Displ , 
(in) 



Y-Displ . 
(in) 



Rotation 

(rad) 



1 


1 


0.00000 


0.00000 


0.00000 


2 


1 


- .04866 


-.16462 


- .00183 


3 


1 


- .01909 


- . 51629 


- .00320 


4 


1 


.05002 


-.99971 


- . 00413 


5 


1 


. 12996 


- 1.56879 


- .00470 


6 


1 


.19793 


- 2.17521 


-.00482 


7 


1 


. 22448 


- 2.48533 


- .00506 


8 


1 


.24696 


- 2.80579 


-.00520 


9 


1 


. 26556 


- 3 . 41666 


-.00460 


10 


1 


.24262 


- 3.93464 


- .00372 


11 


1 


.23091 


- 4.04874 


. -.00366 


12 


1 


.20001 


- 4.28527 


- .00331 


13 


1 


. 18682 


- 4.36562 


- .00322 


14 


1 


. 10299 


- 4.66788 


-.00171 


15 


1 


0.00000 


- 4.77585 


0.00000 


16 


1 


-.10299 


- 4.66788 


.00171 


17 


1 


- .18682 


- 4.36562 


.00322 


18 


1 


-.20001 


- 4.28527 


.00331 


19 


1 


-.23091 


- 4.04874 


.00366 


20 


1 


-.24262 


- 3.93464 


.00372 


21 


1 


- . 26556 


- 3.41666 


.00460 


22 


1 


-.24696 


- 2.80579 


.00520 


23 


1 


-.22448 


- 2.48533 


.00506 


24 


1 


- . 19793 


- 2.17521 


.00482 


25 


1 


-.12996 


- 1.56879 


.00470 
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CRITICAL LOAD OM PARABOLIC ARCH 



WHITNEY'S ARCH 



lase Results 
loint Load 


X-Displ. 


Y-Displ. 


Rotation 


dumber 


Case 


( in) 


( in ) 


(r.at5) 


26 


1 


- . 05002 


- .99971 


.00413 


27 


1 


. 01909 


- . 51629 


. 00320 


28 


1 


. 04866 


- . 16462 


. 00183 


29 


1 


0.00000 


0.00000 


0.00000 



8G 



lase Results 
oint. Load 

fvunber Case 

1 1 

2 1 

3 1 



I* 1 

5 1 



6 1 



7 1 



8 1 

9 1 



10 1 



11 1 

12 1 

13 1 

14 1 

15 1 

16 1 

17 1 

18 1 

19 1 

20 1 

21 1 

22 1 

23 1 

24 1 

25 1 



CRITICAL LOAD ON WHITNEY ARCH 
*♦♦ JOINT DISPLACEMENTS *** 



X-Displ. Y-Displ. 



( in) 


(in) 


0.00000 


0.00000 


-.33012 


.40387 


- .89038 


1.37969 


-1.45051 


2.48996 


-1.87802 


3.39780 


-2.12684 


3.87432 


-2.18403 


3.89519 


-2.19963 


3.74565 


-2.14320 


2.96761 


-2.03244 


1.59980 


-2.00495 


1.16961 


-1.95114 


.10760 


-1.93647 


-.31791 


-1.91131 


-2.61165 


-1.99977 


-5.01194 


-2.18664 


-7.23630 


-2.43349 


-9.00986 


-2.48345 


-9.29158 


-2.61526 


-9.92319 


-2.66847 


-10.13806 


-2.82642 


-10.50458 


-2.83337 


-10.06358 


-2.76211 


-9.56067 


-2,63772 


-8.88741 


-2.21944 


-7.08717 


87 





Rotation 

(rad) 

0.00000 

.00688 

.00962 

.00913 

.00628 

.00194 

- .00086 

- .00373 

- .00885 

- .01357 

- .01479 

-.01710 

- .01787 

-.01989 

- .01969 

-.01708 

-.01235 

-.01129 

-.00805 

-.00654 

.00018 

.00679 

.00969 

.01237 

.01691 



CRITICAL LOAD ON WHITNEY ARCH 



’ase Results 








foint 


Load 


X-Displ. 


Y-Displ . 


Rotation 


lumber 


Case 


( in) 


( in ) 


(rad) 


26 


1 


-1.59238 


-4.85981 


.01892 


27 


1 


-.86867 


-2.61023 


.01736 


1 28 


1 


-.23086 


-.77844 


.01128 


29 


1 


0.00000 


0.00000 


0.00000 



CRITICAL LOAD ON WHITNEY ARCH 



JOINT DATA ♦** 



X - coord. 


Y - coord. 


X - Degree 


Y - Degree 


Z - Degree 


( feet) 


( feet ) 


of Freedom 


of Freedom 


of Freedom 


0 


0 


0 


0 


0 


9.919 


4.963 


1 


1 


1 


19.87 


9.438 


1 


1 


1 


29.84 


13.43 


1 


1 


1 


39.83 


16.92 


1 


1 


1 


49.84 


19.92 


1 


1 


1 


54.86 


21.24 


1 


1 


1 


59.88 


22.42 


1 


1 


1 


69.92 


24.43 


1 


1 


1 


79.94 


25.95 


1 


1 


1 


82.45 


26.25 


1 


1 


1 


87.96 


26.8 


1 


1 


1 


89.98 


26.96 


1 


1 


1 


100 


27.48 


1 


1 


1 


110 


27.5 


1 


1 


1 


120 


27.08 


1 


1 


1 


129.98 


26.22 


1 


1 


1 


131.98 


26 


1 


1 


1 


137.46 


25.31 


1 


1 


1 


139.95 


24.97 


1 


1 


1 


149.94 


23.31 


1 


1 


1 


159.92 


21.22 


1 


1 


1 


164.88 


20.02 


1 


1 


1 


169.86 


18.72 


1 


1 


1 


179.83 


15.82 


1 


1 


1 


189.89 


12.49 


1 


1 


1 


199.87 


8.774 


1 


1 


1 


209.99 


4.583 


1 


1 


1 


220 


0 


0 


0 


0 



Degree of Freedom: 0=restrained l=free i=coupled to joint 'j' 



y(s: I 




Kelley 

Safety and stability 
concrete barrel shell 
Structures. 



in 

roof 



L 



Thesis 

K264 Kelley 

c.l Safety and stability in 

concrete barrel shell roof 



structures . 



